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We investigate transport properties of quantized chaotic systems in the short wavelength limit. We 
focus on non-coherent quantities such as the Drude conductance, its sample-to-sample fluctuations, 
shot-noise and the transmission spectrum, as well as coherent effects such as weak localization. We 
show how these properties are influenced by the emergence of the Ehrenfest time scale te. Expressed 
in an optimal phase-space basis, the scattering matrix acquires a block-diagonal form as te increases, 
reflecting the splitting of the system into two cavities in parallel, a classical deterministic cavity 
(with all transmission eigenvalues either or 1) and a quantum mechanical stochastic cavity. This 
results in the suppression of the Fano factor for shot-noise and the deviation of sample-to-sample 
conductance fluctuations from their universal value. We further present a semiclassical theory for 
weak localization which captures non-ergodic phase-space structures and preserves the unitarity of 
the theory. Contrarily to our previous claim [Phys. Rev. Lett. 94, 116801 (2005)], we find that 
the leading off-diagonal contribution to the conductance leads to the exponential suppression of the 
coherent backscattering peak and of weak localization at finite te. This latter finding is substantiated 
by numerical magnetoconductance calculations. 

PACS numbers: 73.23.-b, 74.40.+k, 05.45.Mt, 05.45.Pq 



I. INTRODUCTION 

Closed chaotic systems are classically characterized by 
ergodicity, mixing and a positive Kolmogorov-Sinai (KS) 
entropy 1] . These three characteristics form a hierarchy: 
mixing systems are ergodic, and systems with positive 
KS entropy are mixing, but the reverse is not necessarily 
true. Ergodicity means that phase-space averages equal 
time averages, while the definition of both mixing and 
KS entropy requires the introduction of some phase-space 
coarse-graining. For mixing, one needs to define finite- 
sized phase-space cells inside which points originating 
from two initially well separated distributions of initial 
conditions are equally likely to be found. As time goes 
by, mixing occurs on smaller and smaller scales, i.e. the 
minimal volume of these cells decreases. The KS entropy 
is defined from the measure of the intersection of the cells 
with their back evolution under the system dynamics. A 
positive KS entropy means an exponential production of 
information, and thus the generation of randomness in 
the Kolmogorov sense, as more and more different tra- 
jectories emerge from apparently indistinguishable initial 
conditions 0. For closed systems, the KS entropy is re- 
lated to the exponential sensitivity to initial conditions, 
and equals the sum of the associated positive Lyapunov 
exponents 0,0- 

The situation becomes different once the system is 
open and scattering trajectories are considered. Phase- 
space structures emerge then which are prohibited by 
ergodicity and mixing, even in systems which have a pos- 
itive KS entropy when closed. These structures and their 
influence on quantum transport are the focus of this ar- 
ticle. We will see how their occurrence affects transport 



through open quantized chaotic systems in the semiclas- 
sical, short wavelength limit. They result in a splitting of 
the cavity into two cavities in parallel, one where trans- 
port is ruled by classical determinism and one where 
transport exhibits quantum stochasticity. 



A. Classical chaos in open systems 

We specialize to two-dimensional chaotic cavities in a 
two-terminal geometry. Typical nonergodic structures 
occurring in such open chaotic systems are illustrated 
in Fig. ^ A color plot of the transmission probability 
is shown on a phase-space projection of one of the two 
openings. The horizontal axis gives the position on a 
cross-section of the opening, normalized by the cavity 
perimeter L, and the vertical axis gives the momentum 
component of injection into the system, parallel to the 
cavity boundary, and normalized by the Fermi momen- 
tum py- Both the real-space set-up and the dimensionless 
phase-space we use are defined in Fig. [5] 

Band-like structures such as those appearing in Fi&JD 
have been reported and discussed earlier 0, 0, IE IE 00] 
All scattering trajectories whose initial point lies in one 
of the bands have approximately the same dwell time 
through the system The typical dwell time r of a 
band determines its width as ~ (W/L) exp[— Ar] £| (W 
is the width of the opening and A is the Lyapunov expo- 
nent). Thus the largest blue and red bands in Fig. 2] 
respectively correspond to direct reflection and trans- 
mission, while thinner bands correspond to longer dwell 
times through the system. Trajectories inside a band are 
transported in one bunch, and the phase-space volume 
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FIG. 1: (Color online) Classical phase-space color plot of 
the transmission probability from the phase-space projection 
of the injection lead (see Fig. The phase-space has been 
coarse-grained by a rectangular grid with 9 ■ 10 4 cells, and the 
transmission probability in each cell has been calculated by 
time-evolving 10 4 classical trajectories per cell. Higher trans- 
missions correspond to red, lower transmissions to blue. The 
dynamical system used is the open kicked rotator as defined 
in Section fyTl 



FIG. 2: (Color online) Sketch of a two-terminal open chaotic 
cavity (top left) and its phase-space represented as a Poincare 
surface of section on the boundary of the cavity (bottom 
right). All momenta on the energy surface (E = p F /2m) 
are parametrized by the tangential momentum p running 
from — pf to pf. All possible positions on the boundary are 
parametrized by q running from zero to L, where L is the 
circumference of the cavity. The phase-space is made dimen- 
sionless by normalizing momenta and real-space coordinates 
with pf and L respectively. 



B. Quantum chaos in open systems 



they occupy is blocked for other trajectories by Liou- 
ville's theorem. Because trajectories remain inside the 
system for a finite time, the definition of ergodicity, that 



tt 1 / dpdq /(p,q;t) = lim r 1 / dt/(p , q ; t), 
Jn T ^°° Jo 

(1) 

for almost all functions /(p,q;i) and almost all phase- 
space points (po,qo) G f2 no longer holds, but depends 
on (po,qo)- The time-integral on the right-hand side of 
Eq. extends only up to the dwell time to of the one tra- 
jectory going through (po,qo), and accordingly Eq. 
cannot be preserved over the full phase-space Q. Simul- 
taneously, mixing occurs on a given scale only for the 
subset of trajectories longer than some finite dwell time. 
Scattering trajectories through open systems have a con- 
tinuous distribution of dwell times P(t) and because of 
the exponentially decreasing volume of scattering bands, 
mixing occurs on exponentially smaller scales on longer 
trajectories. Reversing the argument, a given phase- 
space resolution volume £ corresponds to a time scale 
ps A -1 ln[(VK/L) 2 /^]. Long trajectories with t > ex- 
hibit mixing on the scale £, while short trajectories with 
t < lie on bands well resolved by £-cells, which there- 
fore do not have the mixing property. These two sets of 
classical scattering trajectories have no phase-space over- 
lap. 



A finite resolution scale emerges naturally when the 
system is quantized: the phase-space becomes tiled with 
cells of volume 2irh. For particles with a finite Fermi 
wavelength Ap, this is equivalent to tiling the dimen- 
sionless phase-space projection of the leads (see Fig. |2J) 
with cells of volume h c g = Af/L, the effective Planck's 
constant. This leads to the existence of a finite num- 
ber N = IntpVF/Ap] of conduction channels through the 
system. As h c s is made smaller and smaller, all classical 
parameters being kept constant (the semiclassical limit), 
more and more of the band structures of Fig. ^ are re- 
solved (see Fig. [SJ. Conversely, more and more of the 
conduction channels are supported by one and only one 
classical transmission or reflection band, and thereby be- 
come deterministic. It is thus natural to investigate the 
effect that the lack of mixing of short trajectories has on 
properties of open quantum chaotic systems. 

Transport through ballistic quantum cavities, so-called 
quantum dots, has been investigated intensively in recent 
years 10]. In the regime where the dot's size is much 
larger than the Fermi wavelength, L > Af, transport 
has been shown to depend on the integrability or lack 
thereof of the classica l dy namics, as determined by the 
confinement potential |l2j . Most experimental inves- 
tigations so far have focused on the limit of few, < 10 
conduction channels, where it has been found that quan- 
tum transport in the chaotic case exhibits a universality 
which is well captured b y th e Random Matrix Theory 
(RMT) of transport [Hill. Recently, a semiclassical 
theory for the conductance of ballistic cavities has been 
developed , confirming the common belief that RMT 
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universality applies at least to a certain regime of ballistic 
chaos. 

It is however well-known that physically relevant time 
scales restrict the range of validity of RMT. In recent 
years it has become clear that the Ehrenfest time te 
does this in ballistic quantum chaotic systems, with RMT 
ceasing to be valid when te becomes relevant. The 
Ehrenfest time is the time it takes for the chaotic clas- 
sical dynamics to stretch a narrow wave packet, of spa- 
tial extension Af, to some relevant classical length scale 
C. Since the stretching is exponential in a chaotic sys- 
tem, one has te = A -1 \n[£/h c ffL] 16(. The scattering of 
an initially localized wavepacket into all possible modes 
(similar to s-wave scattering on a restricted portion of 
phase-space) is only established after classical mixing has 
set in on the scale /i e ff, i.e. for times longer than te- 
For shorter times, the quantum dynamics is determinis- 
tic. One thus expects deviations from RMT to emerge 
as te/tq increases. 

This line of reasoning has been qualitatively confirmed 
in the cavity transport experiments of Ref . [XtJ , which ob- 
served a significant reduction of the electronic shot-noise 
power below its RMT value upon opening the cavity more 
and more. This reduction is very likely due to an increas- 
ing fraction of deterministic channels in the transmission 
spectrum of the cavity which can be understood as fol- 
lows. The shot-noise power (the intrinsically quantum 
part of the fluctuations of the electronic current) is usu- 
ally quantified by the dimensionless Fano factor F, the 
ratio of the shot-noise to the Poissonian noise [18( , which 
can be expressed in term of the transmission spectrum 
Ti € [0, 1] of the cavity as [H| 

^ = 5>(l-T0/5>. (2) 

i i 

Hence deterministic channels, those having Ti = or 1, 
do not contribute to F. Such channels appear as the 
classical bands discussed above are quantum mechani- 
cally resolved, which can be achieved either by reduc- 
ing the dwell time or by reducing the Fermi wavelength. 
The former change was more appropriate for the experi- 
mental set-up of Ref. 01 i an( A the latter change is illus- 
trated in the numerics of Fig. |3] We show three quantum 
phase-space plots for a fixed classical set-up (the same as 
for Fig.^). Plotted is the quantum transmission proba- 
bility ((x,p)\T\(x,p)} for a fixed grid of initial coherent 
states \(x,p)), i.e. isotropic Gaussian wavepackets cen- 
tered on (x,p). The three panels from left to right cor- 
respond to smaller and smaller Ap. It is seen that as Ap 
decreases, finer and finer phase-space structures are re- 
solved. Moreover, coherent states entirely lying on deep 
red (deep blue) regions have ((x,p)\T\(x,p)) = 1 (0), and 
are therefore eigenstates of T. Any of them can carry a 
quantum channel which does not contribute to shot-noise 
(the total number of deterministic channels is obtained 
only after the orthogonalization of the coherent states, 
see below). With decreasing Af, the number of deter- 
ministically transmitted coherent states increases faster 



than the total number of channels, inducing the reduc- 
tion of the shot-noise power below its RMT value. 

The suppression of the Fano factor in the semiclas- 
sical limit was anticipated long ago 0|. More re- 
cent works quantitatively predicted a suppression F cx 
exp(— te/td), in term of the new Ehrenfest time scale and 
the average dwell time td through the system [20l l2l| , 
a suppression which was related to the phase-space reso- 
lution picture of Ref. k| given above and confirmed nu- 
merically in Refs. Ref- H3 presented a phase- 
space semiclassical approach resolving the classical bands 
which showed that the fraction of deterministic trans- 
mission eigenvalues not contributing to noise is cx [1 — 
exp(— te/td)]. Following the numerous recent investiga- 
tions of the quantum-to-classical correspondence in open 
systems, which we now proceed to briefly summarize, it 
has become clear that te/td — * is a necessary condition 
for complete RMT universality [23|. As is illustrated in 
Fig- El this condition is never satisfied in the semiclassical 
limit H c ff — > 0. 

Following Ref. [24| which suggested that the existence 
of a finite te discriminates quantum chaotic from quan- 
tum disordered systems (with the former class exhibiting 
deviations from universality), many investigations have 
been devoted to the study of open quantum chaotic sys- 
tems at finite te |23|. Focusing on transport, it is by now 
well established numerically that, as te/td — * oo, the 
Fano factor disappears [j| , and sample-to-sample conduc- 
tance fluctuations lose their universality. Simultaneously, 
parametric conductance fluctuations appear to remain 
universal 0, 0] ■ The weak localization correction to the 
conductance was first predicted to disappear exponen- 
tially as Sg oc exp[— 2te/td] jlij or Sg oc exp[— te/td] 
25]. Early numerical works, on the other hand, con- 
cluded that Sg is independent of te/td [2(|. This puz- 
zle was resolved in Ref. 27] (though some skepticism re- 
mained, on the part of the authors of the present arti- 
cle amongst others) whose analytical approach predicted 
Sg oc exp[— te/td] in agreement with numerical simula- 
tions. Below we will present both an analytical, semi- 
classical treatment of weak localization with a special 
emphasis on backscattering and unitarity, and numerical 
magnetoconductance data giving a microscopic confirma- 
tion of the conclusion Sg oc exp[— te/td] of Ref. |27| . 

There are currently several theories for open quantum 
chaotic systems at finite te- First, the stochastic qua- 
siclassical theory mimics the post-Ehrenfest time mode- 
mixing by introducing fictitious random scatterers with 
a scattering rate appropriately tuned [23, l24l |27| . It 
is developed from standard methods in disordered sys- 
tems, but breaks time-reversal symmetry at the classi- 
cal level already. Second, there is a semiclassical theory 
|25| , which until now had not been shown to preserve the 
unitarity of the scattering matrix (and hence conserve 
the current). Third, a phenomenological model origi- 
nating from Ref. models the total electronic fluid as 
two coexisting phases, a classical and a quantum one. 
At this level, the theory is referred to as the two-phase 
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FIG. 3: (Color online) Quantum phase-space color plots of the transmission probability from the injection lead. The system 
is the quantum equivalent of the classical system of Fig. Q The phase-space has been coarse-grained by a rectangular grid 
with 9 • 10 4 points. Starting from each point of the grid, an isotropic Gaussian wave-packet has been time-evolved and its 
transmission probability calculated. From left to right, the three panels correspond to decreasing effective Planck's constant 
h cS = 2-r/M with M = 512 (with a conductance g = 22.4 and a Fano factor F = 0.193), 8192 (g = 375.9 and F = 0.121) and 
131072 (g — 5990.8 and F ~ 0.08) respectively. More and more fine-structured details of the classical phase-space are resolved 
as h e g — > 0. Higher transmissions correspond to red, lower transmissions to blue color. 



Two-cavity regime 




FIG. 4: (Color online) Schematic of the different transport 
regimes through a ballistic chaotic cavity with perimeter L 
coupled to leads of width W . Above the separatrix, the sys- 
tem splits into two effective cavities, one purely classical and 
the other quantum mechanical. Only the latter contains quan- 
tum effects such as non-deterministic transmission and quan- 
tum interferences. Below the separatrix W = (ApL) 1 ' 2 (solid 
curve), all modes are quantum mechanical (grey area). The 
red dashed arrows indicate the semiclassical limit of h e « — > 
at fixed classical parameters. The slope of the arrows is given 
by the inverse dwell time Tb. 



fluid model 0, • With the additional surmise made 
in Ref. Q that the quantum phase has RMT properties, 
one gets the effective RMT model. The effective RMT 
model successfully explains the behavior of shot-noise, 
the transmission spectrum and conductance fluctuations, 
but is in contradiction with the disappearance of weak lo- 
calization 27] . The suppression of weak localization at 
large te/td invalidates the effective RMT model, but not 
the two-phase fluid model, as we will explicitly see below. 



C. Outline of this article 

In this article, we focus on quantities such as the aver- 
age conductance, shot-noise and the transmission spec- 
trum in ballistic chaotic cavities, as well as the weak 
localization correction to the conductance as H e g — > 
0. These quantities are strongly influenced by the 
emergence of the open cavity Ehrenfest time Tg P = 
A -1 \n[(W/L) 2 /2Trh e s] (28|. All classical parameters be- 
ing fixed, that limit inevitably turns any system into a 
nonuniversal quantum chaotic one as t^P becomes finite 
(see Fig.0J. We calculate the scattering matrix in a basis 
that optimally resolves phase-space structures and show 
that the system splits into two cavities in parallel. This 
provides a foundation for the two-phase fluid model. We 
go significantly beyond our previous work, Ref. , with 
(i) a detailed construction of a basis which optimally re- 
solves those phase-space structures and (ii) an explicit 
semiclassical calculation of the weak localization correc- 
tion to the conductance which preserves the unitarity of 
the scattering matrix to leading order. 

The outline of the paper goes as follows. In Section HTl 
we discuss the nonergodic classical structures of open 
chaotic systems such as those shown in Fig. ^ and the 
Ehrenfest time scale that accompanies them. Our task 
requires that we resolve quantum mechanically these clas- 
sical phase-space structures. This suggests that we em- 
ploy a semiclassical theory. The existing such theories 

mm 

have to be replaced by a phase-space resolving 
theory, which requires the construction of an appropri- 
ate orthogonal mode basis. This is done is Section ITTT1 
We then write the system's scattering matrix in this ba- 
sis, and show how this results in phase-space splitting at 
the quantum level in Section IIVI Amplifying on that, 
we show how deterministic modes emerge and calculate 
the average Drude conductance and its sample-to-sample 
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fluctuations at large t^ p . In Section we present our 
semiclassical theory for the weak localization correction 
to the conductance. Several key aspects absent in previ- 
ous semiclassical treatments are stressed here. In partic- 
ular we present the first semiclassical calculation of the 
coherent backscattering peak at finite t^/td- We show 
that both coherent backscattering and weak localization 
are exponentially suppressed oc exp[— t^/to] with the 
closed cavity Ehrenfest time r^ 1 = A -1 ln[?Cf ]. This is 
so, because weak localization and coherent backscatter- 
ing come from trajectories longer than Tg P + Tg , which 
have an exponentially small relative weight exp[— /tq] 
in the stochastic block of the scattering matrix. The 
existence of two separated fluids is however confirmed. 
We demonstrate that the classical phase-space structures 
(which give rise to phase-space splitting) must be in- 
cluded in the semiclassical treatment to preserve the uni- 
tarity of the scattering matrix at nonvanishing value of 
t"e P /td- We also point out the origin of the discrepancy 
between our final conclusion, that 5g oc exp[— t^/td], 
and our earlier claim that Sg remains at its universal 
value in the deep semiclassical limit j^]. We finally 
present numerical magnetoconductance data confirming 
our theory. Summary and conclusions are presented in 
Section IVlTl and technical details discussed in the Ap- 
pendices. 



II. CLASSICAL SCATTERING STRUCTURES 
AND EHRENFEST TIMES 

A. Transmission and reflection bands 

We consider classical scattering trajectories. They are 
injected into the cavity from one of the two leads, say 
the left (L) lead, with initial position q and momentum 
p on a cross-section of the lead. The momentum is on 
the Fermi energy surface E = pp/2m. The trajectory 
is determined by ballistic motion inside the confinement 
potential defining the cavity, until the particle hits the 
boundary between the cavity and one of the leads, at 
what time it escapes. Throughout this paper we will 
consider perfectly transparent leads. 

Scattering trajectories are not isolated, instead they 
occur in bands (see Fig. ^ . As mentioned in the intro- 
duction, a scattering band is a phase-space structure oc- 
curring in open systems, even when their closed counter- 
part is fully chaotic. It contains a set of trajectories which 
exit at approximately the same time through the same 
lead 0, having followed similar trajectories through the 
cavity, in the sense that any trajectory in the band can 
be topologically deformed into any other. The situation 
is sketched in Fig. [5] Bands on the injection lead are de- 
fined by the overlap of that lead with the time-reversed 
evolution of the exit lead, including absorption at both 
leads. 

For an individual system, the exact number and area of 
the bands depends on the specifics of the lead positions 
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FIG. 5: (Color online) Time evolution of an L to R trans- 
mission band (parallelograms) on the boundary of the cavity 
of Fig. |21 The leads are indicated by the shaded rectangles. 



and widths, and the cavity shape. However, averaged 
properties calculated over an ensemble of cavities with 
the same dwell time and Lyapunov exponent can be cal- 
culated. The asymptotic average survival probability is 
exponential |29j . 



p(r) = exp[-r/r D ]- 



(3) 



It depends solely on the average dwell time td given by 



ttA 



v F (W L + W R ) 



(4) 



Here A is the area of the cavity and we considered leads 
with different widths Wl.r- 

In the dimcnsionless phase-space defined in Fig. |3 
where momenta and distances are measured in units of 
Pf and L respectively, the injection lead (we always as- 
sume this is the L lead) has a dimensionless phase-space 
area of 



2WL/L. 



(5) 



The fraction of the lead phase-space which couples to 
transmitting trajectories is ~ Wr/(Wl + Wr), with 
the rest coupling to reflecting trajectories. The average 
phase-space area of a single transmission (L— >R) or re- 
flection (L— >L) band which exits at time t, is given by 



(r)> 



W L W K 



-Xt 



(6) 



where K = L,R. The average number (riL— >k(t)) 
of bands exiting at time r through lead K is 
given by multiplying El/(£l— >k(t)) by the probability 
(Wk/L) exp[— t/td] to escape through the K lead at r. 
Hence one has 

(n L _ L (r)) = (n L ^ R (r)) ~ exp[Ar(l - (Atd)" 1 )]. (7) 

Since we assume fully developed chaos, Atd 3> 1, we see 
that the average number of bands diverges as r goes to 
infinity even though the sum of their phase-space areas 
goes to zero. We also note that the average number of 
reflection and transmission bands are equal, with only 
their areas being dependent on Wi, r,. 
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B. Ehrenfest times and modes on classical bands 

The Ehrenfest time scale emerges out of the quantum 
mechanical coarse-graining of phase-space and the par- 
tial resolution of scattering bands. It is the time it takes 
for quantum mechanical uncertainties to blow up to some 
relevant classical scale L in chaotic systems. The scale 
C depends on theproblem at hand, e.g. on whether the 
system is closed , or open [23, 123 • For the transport 
set-up we will focus on, this scale is related to the area 
of scattering bands. Large scattering bands, those with 
phase-space area greater than 2irH, can carry a number 
of modes of order their phase-space area divided by 2-Kb,. 
All those modes are classical, deterministic and exhibit 
no quantum effects. They are supported by trajecto- 
ries shorter than the Ehrenfest time. The small bands 
on the other hand, those with area less than 2irh carry 
less than a full mode, which generates quantum (stochas- 
tic) modes, sitting on many small bands with dwell times 
longer than the Ehrenfest time, and hence being partially 
transmitted and partially reflected. Eq. © then defines 
two open cavity Ehrenfest times for states injected from 
the L lead, one for transmitting bands and one for re- 
flecting bands 



% K = A" 1 In 



W L W K 



2Tth c ffL 2 



K = L,R. 



(8) 



The difference between and is only logarith- 
mic in Wr/Wl- We will often neglect it and con- 
sider instead the symmetric open cavity Ehrenfest time 
t e = A" 1 \n[(W/L) 2 /2nh cS }. The open cavity Ehren- 
fest times, t% can be interpreted as the time it takes 
for a wavepacket of width Wl, /L along the stable mani- 
fold of the hyperbolic classical dynamics to evolve into a 
wavepacket with width Wk/L in the unstable direction. 

We can readily estimate the number of quantum scat- 
tering modes. The proportion of the L lead phase-space 
which couples to trajectories to the K lead with r > t^ k 

is on average e~ T ^ K / Td Wk / (Wl+Wr). Thus the average 
number of quantum modes in the L lead is 

LL / LR / 

/v T e~ r E /to 1 /Vne~ TE ' TD 

w>-*- — ■ m 

All other modes of the L lead are in transmission bands 
with r < Te R or reflection bands with r < Tg L , and so 
they are all classical modes. Their number is thus 



= N L (l 



0[(Xtd) 



(10) 



The counting argument leading to these estimates finds 
a rigorous derivation below in section IIII Bl where we 
explicitly cover scattering bands with an orthonormal 
phase-space (PS) basis. There, we also identify a third 
class of states, overlapping significantly but still only par- 
tially with large bands with r < Tg K . Because of their 
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FIG. 6: (Color online) Main plot: Plot of the real-space 
wavefunction components |(Q|ps; 0, 0)| 2 of a PS-basis state 
(black line), and |(Q|cs; 0, 0)| 2 of a coherent state (red line). 
Both wavefunctions are symmetric in the dimensionless phase- 
space under Q <-> P, up to a scaling factor. Inset: Logarith- 
mic plot of |(Q|ps; 0, 0)| 2 , showing that the PS-state decays 
exponentially in position (dashed line). The oscillations in 
I (Q|ps; 0, 0}| 2 ensures it is orthogonal to PS-states centered 
at finite Q, while its broader central peak ensures that it is 
orthogonal to PS-states centered at finite P. 



small relative number however, these states only have a 
subdominant effect on the system's properties. 

Note that the number of classical PS-states goes like 



^eff^ ATD ^)' wnne t ne number of quantum PS- 



-1/(Atd) 



In the semiclassical limit 



states goes like h~ s 
^■off — ► we see that the number of quantum PS-states 
goes to infinity, while their fraction goes to zero. 

III. THE PHASE-SPACE BASIS 

In order to formally split classical from quantum 
modes, our task now is to construct a complete orthonor- 
mal basis resolving maximally the scattering band struc- 
ture of the classical phase-space. This requires an op- 
timal resolution in both real-space and momentum co- 
ordinates. To achieve that, we use results from wavelet 
analysis. 

A. Existence of orthogonal phase-space bases 

The existence of complete orthonormal bases with 
states exponentially localized in time and frequency has 
been proven in the context of wavelet analysis [3l|. We 
use such a basis as a PS-basis, in which each basis state is 
exponentially localized in position and momentum. We 
are unaware of any such basis which has closed form ex- 
pressions for the basis states. There are however numer- 
ous algorithms for generating such bases [3l|]. In Ap- 
pendix [B] we give such an algorithm which iteratively 
orthogonalizes a complete but non-orthogonal basis of 
coherent states, generating a set of PS-states of the form 



7 



shown in Fig. EJ While we give this explicit example, we 
emphasize that our theory only requires the existence of 
such a basis. We use the fact that each basis state is 
exponentially localized in position and momentum, and 
that any such complete orthonormal basis remains com- 
plete and orthonormal under any rotation, translation 
or area-preserving stretch in phase-space. Having con- 
structed the PS-basis, the transformation which takes us 
from lead modes to PS-states is unitary since both bases 
are orthonormal. 



B. The optimal phase-space basis 

In a recent letter j^, we constructed an orthonormal 
PS-basis on a square von Neumann lattice. This basis 
is simple to explain and work with, however, it underes- 
timates the number of classical modes and in particular 
leads one to predict that the open cavity Ehrenfest time 
is half its correct value 32] . To obtain the correct value 
of te, the von Neumann lattice must be adapted to fit 
in the classical band structures as best as it can. This is 
done band by band. For parallelogram bands, the pro- 
cedure is to rotate and stretch the originally square von 
Neumann lattice to a parallelogram von Neumann lat- 
tice. This is illustrated in Fig. Each lattice cell still 
covers an area 2?:h and the intraband orthogonality is 
ensured. The interband basis orthogonality is preserved 
due to the exponentially small overlap of PS-basis states 
in different bands (classical bands do not overlap thanks 
to Liouville's theorem; this effect has been termed Liou- 
ville blocking in Ref. j22|), except for a minority of states 
lying directly at the boundary of the band which we will 
deal with below. This procedure can still be applied as 
long as the band's curvature is not too large, or for bands 
which look more trapezoidal than parallelogram-like. All 
one needs to do then is adapt locally the aspect ratio of 
the von Neumann lattice, as shown on Fig. [8] Bands with 
small curvatures dominate at short dwell times. How- 
ever, some bands with larger dwell times inevitably dis- 
play a fold. For those bands, the procedure is to bend 
the von Neumann lattice along the axes defined by the 
eigenvectors of the stability matrix of the classical dy- 
namics at each point in the band's phase-space, as in 
Fig. |S| The aspect ratio of the lattice is chosen so it 
obeys Eq. I|llf) locally. For intermediate values of h e g, 
the local curvature of the resulting lattice destroys the 
orthogonality of the PS-states, however, as H e s is reduced 
further and further, the curvature drops out of the prob- 
lem, any smooth curve being locally well approximated 
by a straight line. To formally show that an optimal 
orthogonal basis can be generated from the square von 
Neumann basis of Ref. [22| it is thus sufficient to (i) con- 
sider parallelogram and trapezoidal bands only, keeping 
in mind however that (ii) at any finite value of ft e ff, de- 
viations from parallelogram shape generates additional 
non-optimal PS-basis states. The latter result from a 
further orthogonalization (e.g. Gram-Schmidt) required 




FIG. 7: (Color online) Sketch of two parallelogram bands 
(green areas) with the PS-states superimposed on them (el- 
lipses). The upper band has a phase-space area which is a 
bit more than sixteen times 2irh, while the lower band has a 
phase-space area a bit more than 2irh. In both cases the basis 
is optimized, the lattice of PS-states is stretched and rotated 
such that the maximum number of PS-states can be htted 
into each band (solid-edged ellipses) with the minimum num- 
ber partially in the band (dashed-edged ellipses). Note that 
these optimally chosen PS-states have the same aspect ratio 
as the classical band in which they sit, thus their longitudinal 
and transversal extensions, x and y, are given by Eq. JTT1 . 



for edge-of-band states and for states on folds whose cur- 
vature is not yet well resolved at this value of h c g . Below 
in Appendix El we will see that those states build up 
a negligible fraction iV qm /(ArD) <C iV qm <C N c \ of the 
total number of modes in the semiclassical limit. The 
completeness of the basis follows from the orthogonality 
and the conservation of the total number of basis states, 
the above procedure being area-conserving. We are now 
ready to extend the discussion of Ref. [z| and derive an 
optimal phase-space basis for parallelogram bands, which 
gives the correct Ehrenfest time. This is done in a four 
step process. 

Step [i] . Pick a scattering band with phase-space area 
larger than 2irh, and cover it with a lattice of PS-states. 
Both the lattice and the states must be stretched and 
rotated to the same angle and aspect ratio as the band, 
and positioned in such a way as to minimize the number 
of edge-of-band states. This is illustrated in Fig. [7] This 
can be done without relaxing either the mutual orthog- 
onality, or the normalization of the PS-states. Since the 
PS-states have the same aspect ratio as the band, their 
longitudinal and transversal extensions x and y, as indi- 
cated in Fig. [7| are given by (r is the dwell time of the 
band under consideration) |33| 

(2^ cff WVIEr) 1/2 exp [At/2] ; V 2nh cS /x. (11) 

While we pay attention to minimizing their number, we 
do not include edge-of-band states in the basis at this 
stage. We will deal with them later in step [iv]. 
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FIG. 8: (Color online) Sketch of a trapezoidal (top green 
area) and a folded (bottom green area) scattering band with 
a lattice of PS-states superimposed on them. Both bands 
cover a phase-space area larger than 2irh. In both cases the 
basis is optimized, the lattice of PS-states is locally stretched 
and rotated such that the maximum number of PS-states can 
be fitted into each band, some of them are indicated by solid- 
edged ellipses. 



Step [iv]. Steps [i-iii] generate an incomplete or- 
thonormal basis in the vector-space of lead modes. This 
basis can be made complete by adding the adequate num- 
ber of states orthogonal to those generated in steps [i-iii] . 
This construction gives us very little information about 
the nature of these additional states, except that they 
must sit on more than one band. They thus evolve in 
a quantum, stochastic manner, and we refer to them as 
quantum PS-states. The set of quantum PS-states di- 
vides into two broad categories: firstly those which sit 
on many classical bands, secondly those which sit mostly, 
but not completely on a single band. We already men- 
tioned the second category of edge-of-band PS-states (the 
dashed ellipses in Fig. 0) . In Appendix El we estimate 
that the number of edge-of-band PS-state is 

N cob ~ (Atd)- 1 ^ (12) 

Hence they form a small fraction of the total number of 
quantum modes, and we do not consider them separately 
from pure quantum modes. 

IV. SCATTERING MATRIX IN THE 
PHASE-SPACE BASIS 



Step [ii]. We evolve the states generated in step [i] 
under the cavity's dynamics. The lattice of PS-states 
on an injection band uniquely determines the lattice of 
PS-states on the exit band. All trajectories in the band 
under consideration exit the cavity after a time shorter 
than the open cavity Ehrenfest time. On this time scale 
the quantum dynamics of the PS-states are well approx- 
imated by the Liouville flow [24j (see also Appendix El . 
This is a well-known property of coherent states that can 
be extended to the exponentially localized PS-states that 
we consider here. Thus a PS-state with initial spread of 
AQ in the unstable direction at r = 0, will evolve into 
a PS-state with spread AQ' ~ AQe AT in the unstable 
direction at time r. The initial spread in the stable di- 
rection is AP ~ h/ AQ, and since the quantum dynamics 
is Liouvillian inside classical bands, it is area preserving, 
i.e. AP' ~ h/AQ'. The PS-states are stretched and 
rotated in the same manner as the exiting band, while 
still forming an orthonormal basis on that band. The 
orientation and stretch of the lattice and PS-states are 
given in linear approximation (which eventually become 
valid as h c s — > 0) by the eigenvalues and eigenvectors of 
the stability matrix of the classical dynamics, and thus 
have the same aspect ratio and angle as the exiting band. 
We choose to use these states to cover that band in the 
phase-space. 

Step [iii] . We simply repeat the process in steps [i- 
ii] for each band with area > 2ttH not yet covered by 
PS-states. The crucial point here is that bands cannot 
overlap; hence if we only place PS-states within classical 
bands, Liouville blocking ensures that PS-states in differ- 
ent bands must be orthogonal with exponential accuracy. 



A. Splitting of the scattering matrix and 
deterministic transmission 

By construction, the PS-basis is chosen so that there is 
a one-to-one correspondence between incoming and out- 
going classical PS-states, given that their quantum dy- 
namics can be approximated by the Liouvillian flow. The 
unitarity of the scattering matrix means that the quan- 
tum PS-states remain orthogonal to the classical ones 
as they evolve inside the cavity. Thus, despite the fact 
that quantum PS-states are not well described by the 
Liouvillian flow, they cannot penetrate the regions of 
phase-space containing bands larger than 27rfi, c ff. In the 
PS-basis, each incoming classical PS-state goes to ex- 
actly one outgoing classical PS-state, while each incom- 
ing quantum PS-state goes to multiple outgoing quantum 
PS-states, but no outgoing classical PS-states. Corre- 
spondingly, the scattering matrix in the PS-basis is of 
the form 

The matrix S d is N cl x N cl w hile the matrix S qm is A^ qm x 
A™ m , with N cl = iV£i + ]Sf^ and 2V qro = iV qm + A^ m . 

The matrix S c i has only one non-zero element in each 
row and each column. The modes on the left and right 
of S can be reordered such that the transmission part t c i 
of S c \ is diagonal with all its non-zero elements in the 
first n elements of its diagonal, where n is the number of 
classical transmission modes. Thus we can write 
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where all non zero-elements of t c i are contained in the 
n x n matrix t c i. Doing the same for t^, r c j and r^, we 
write the classical part of the scattering matrix as 



r c i 
t c i 



t' 






t c l 







fcl 






tcJ 







(15) 



/ 



where t c i and t'. 



are n x n 



matrices, f c i is an (A^ 1 



n) x (Nfj — n) matrix and is an (N^ — n) x (AJj 1 — n) 
The matrix t c i is diagonal with elements given 



matrix. 

by 



(t c l)i 



= e i$i c 



(16) 



The matrix f c i also has exactly one non-zero element in 
each row and each column. Its elements obey 



|(?ci)y| = |(?cl)ji| = 



1 when i reflects to j, 
otherwise. 



(17) 



Anticipating discussions to come, we note that coherent 
backscattering is carried here by matrix elements (?ci)ij 
where the j denotes the time-reversal of i. Since the 
classical probability to go from i to its time-reversal is 
equal to the probability to go from i to itself, the num- 
ber of nonzero such matrix elements is determined by the 
probability to sit on a periodic orbit which does not visit 
the contact to any of the two leads except for one point. 
This probability is zero. The absence of reflection matrix 
elements carrying coherent backscattering in f c i qualita- 
tively explains the exponential suppression of coherent 
backscattering. 

We next calculate the transmission matrix, T = t^t. 
The block diagonal nature of S in the PS-basis given in 
Eq. I|13|l. ensures that T has the same structure in that 
basis, hence 



J-cl 

T, 

' tLt„ 



o 



q m 



cl — "cl" cl ' -"-qm — "qm'qm ■ 

From Eq. I|14ll6fl we get the eigenvalues of T c i, 



T t = 



for 1 < i < 



n. 



for 



n < i < N£ l . 



(18a) 
(18b) 

(19) 



This is what we believe is the first proof of a longstanding 
hypothesis, that in the classical limit the vast majority of 
transmission eigenvalues are zero or one |19) . We know 
that there are A^ 1 such classical modes, with the remain- 
ing modes having a quantum nature, making them un- 
likely to have transmission eigenvalues which are exactly 
zero or one. The block-diagonal structure (|18fl of the 
transmission matrix means that the dimcnsionlcss con- 
ductance, g = J2i Ti and the Fano factor for shot noise 
of Eq. (0), can be written as 



.9 = 

F = 



.9cl + flqm, 

.9cl + 5qm 



(20) 
(21) 



where we have introduced the conductance and Fano 
factor for the two cavities (classic and quantum), 

5cl,qm = SiGcl.qm^i an< ^ ^cl,q m = Eiecl.qm 

ri)]/Ei 6 ci,qm T i]- From E q- 01 we see that 



Fcl = 0. 



(22) 



Anticipating the calculation of the average values (gci.qm) 
(see next subsection), 



F = F qm exp[-r£ R /T D ]. 



(23) 



Since F qra < 1, one see there is an exponential suppres- 
sion of F. 

The PS-basis does not give us much information about 
the quantum PS-states. However since each incoming 
PS-state sits on multiple bands, exiting at different times, 
it must couple to multiple outgoing PS-states. It is ex- 
tremely rare for these outgoing PS-states to be all trans- 
mitting (or all reflecting), and thus we expect that the 
vast majority of their transmission eigenvalues lie be- 
tween zero and one, and thus contribute to shot-noise. 
However the fact that quantum and classical PS-states 
exist in two separate sub-blocks of the scattering and 
transmission matrices, see Eqs. (|13|) and I118[K establishes 
the two-phase fluid model ■ 



B. The average Drude conductance 

From the estimates in Section III Bl the ensemble- 
averaged Drude conductance is the sum of the Drude 
conductances of the quantum and classical cavities, 



(S>D = (ffqm) + (Scl), 



(fqm) 
(ffcl) 



JV L qm JVf jV T ,iV, 



K m + K m Al + a r 

N L N R 

A L A R r , ._ TeL h /t 



T iVR LR / 



(24a) 
(24b) 



1-e" 



A L + A R 

Thus the ensemble averaged Drude conductance is, 

(g)D = a l a r /(a l + a r ). 



(24c) 



(25) 



The splitting of the cavity has little effect on (<7)d, even 
though classical modes and quantum modes do not mix. 
For strong asymmetry there is an additional term of or- 
der Nl(Xtu)^ 2 on the right of Eq. I)24c[l . however our 
calculation is not valid to that order because we ignored 
various order (Atd) -1 terms, such as edge-of-band states, 
in Eq. lfT7TJl. 



C. Sample-to-sample conductance fluctuations 

The precise shape, size and number of the nonergodic 
phase-space structures fluctuates from sample to sample. 
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These fluctuations strongly affect g c \. They are of a clas- 
sical nature, and as such they induce the departure of 
conductance fluctuations from their universal behavior 
0, 0] • Indeed, once ft G ff is small enough that quantum 
mechanics resolves the largest scattering band on aver- 
age, the sample-to-sample conductance fluctuations are 
dominated by the band fluctuations. Since each resolved 
band carries a number oc of channels, one expects 
sample-to-sample conductance fluctuations to exceed the 
universal value in the semiclassical limit, 

<y(9d) « Ks » L ( 26 ) 

The above argument predicts the onset for deviations 
of cr(g) from its universal behavior once the largest 
band is quantum mechanically resolved, i.e. for fi e ff < 
(W/L) 2 exp[— Ar ], where r is the minimal escape time, 
of order the time of flight through the cavity. Both this 
onset and the magnitude (|26[1 of the sample-to-sample 
conductance fluctuations have been observed numerically 
00. 




FIG. 9: (Color online) Sketch of the leading order quantum 
correction to the conductance. Trajectory 7I (solid line) is 
injected at Yq = (yo,#o), crosses itself and escapes at Y = 
(y,0). Its first visit to the crossing (the open dot) occurs at 
Ri = (ri,0i), where ri is the position in the cavity, and 
4>i is the angle of the momentum to a reference axis (not 
shown). Trajectory 72 (dashed line) starts and ends at the 
same positions as 7I, however it avoids the crossing. We 
divide 7I into three segments; legl, loop and leg2. 



V. WEAK LOCALIZATION 

We calculate the leading order quantum correction to 
the Drude conductance. Our treatment applies both to 
the universal (te/td — » 0) and deep semiclassical (finite 
te/td) regimes. We present an explicit treatment of the 
coherent backscattering peak showing that our theory 
preserves the unitarity of the scattering matrix, as well 
as a calculation of the magnetoconductance. Thus far we 
have constructed a special basis (the ps-basis) which is 
aligned along the band structures in the classical phase- 
space. This made it easy to calculate the properties of 
the parts of phase-space covered by bands larger than 
2ttH, allowing us to calculate the deterministic transmis- 
sion eigenvalues T = 0,1. However the complexity of 
the quantum modes of that basis, make it difficult to ex- 
plicitly calculate the transport properties of those modes 
|2lj . We therefore return to the lead- mode basis to cal- 
culate weak localization. Unlike previous works we do 
not neglect the classical bands, however. Indeed, our 
semiclassical approach is able to extract the conductance 
(including weak localization) of both the classical and 
quantum cavities. 

A. Drude conductance 

Semiclassically, the transmission matrix reads 0, , 
tji = -{2^hy 1 ' 2 [dy f dy]T(d^/dyo)y 2 

x (M (yo\i) exp[iS y /h + , (27) 

where \i) is the transverse wavefunction of the ith lead 
mode. This expression sums over all trajectories 7 (with 



classical action S 7 and Maslov index starting at yo on 
a cross-section of the injection (L) lead and ending at y 
on the exit (R) lead. Inserting Eq. I|27|) in the Landauer- 
Buttiker formula for the conductance g = Tr[t't], one 
gets a double sum over trajectories, 7I and 72 and over 
lead modes, \n) and |m). We make the semiclassical ap- 
proximation that J2n(y'\ n }( n \y) — ~~ y) 36J . The 
conductance is then given by a double sum over trajec- 
tories which both go from yo on lead L to y on lead R, 




Here, A~ ( — [dp^/dyo]* ■ Reflection, R = Tr[r^r], is 
given by the same expression, with both yo and y on 
lead L. We are interested in quantities averaged over 
variations in the energy or the cavity shape. For most 
{71,72} the phase of a given contribution, SS/h, will os- 
cillate wildly with these variations, so the contribution 
averages to zero. The most obvious contributions that 
survive averaging are the diagonal ones with 7I = 72. 
These contributions give the Drude conductance H25f) . 

We define P(Y, Y ; t)Sy58St as the product of the ini- 
tial momentum along the injection lead, pfcos#0j an d 
the classical probability to go from an initial position 
and momentum angle Yo = (yo,9o) to within (5y,69) of 
Y = (y, 6) in a time within 5t of t. Then the sum over 
all trajectories 7 from yo to y is 




x P(Y,Y ;t) [••■]y„. (29) 

For an individual system, P has <5-functions for all classi- 
cal trajectories. However averaging over an ensemble of 
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systems or over energy gives a smooth function 



(a) Weak-localization for transmission (conductance) 



(i>(Y,Y ;i)> 



Pf cos 9q cos 6 
2(W L + Wr)t d 



expl-t/rp] . (30) 



This latter expression (|3t)|l is valid as long as no restric- 
tion is imposed on the trajectory inside the cavity. Using 
Eqs. (|29|l and i|30[l to calculate the conductance within 
the diagonal approximation, one recovers the Drude con- 
ductance 



Sdiag = gu 
Pdiag 



_NlNr_ 

N L + N B 

N L + N B 



(31a) 
(31b) 



One also sees that at the level of the diagonal approxima- 
tion, there is unitarity. We stress that, unlike in Ref. |37| . 
we do not include coherent backscattering in the diagonal 
contribution, it is dealt with separately below. 



B. Weak localization for transmission 

A pair of trajectories giving the leading order cor- 
rection to the Drude conductance is shown in Fig. [5] 
The trajectories are paired almost everywhere except in 
the vicinity of an encounter Going through an en- 
counter, one of the trajectories intersects itself, while the 
other one avoids the crossing. Thus, they travel along the 
loop they form in opposite direction. It has been shown 
in Ref. j3g] that for any self-intersecting trajectory with 
a small enough crossing angle e, there exists a partner, 
crossing-avoiding outer trajectory. For the relevant case 
of small e, the probability to find a weak localization pair 
is thus given by the probability to find a self-intersecting 
trajectory. The two trajectories are always close enough 
to each other that their stability is the same, i.e. one can 
set J2-yi 7 2 ^71^72 — * S-yi^i- To evaluate the weak 
localization correction to conductance, we perform a cal- 
culation similar to Ref. [37| , adding the crucial fact that 
pair of trajectories such as depicted in Fig.[5]have highly 
correlated escape probabilities due to the presence of an 
encounter jS^. The situation is depicted in more detail 

in Fig. nni 

The presence of the encounter introduces two new in- 
gredients, both of these were overlooked in Ref. [37| . 
First, pairs of trajectories leaving an encounter escape 
the cavity in either a correlated or an uncorrelated way. 
Uncorrelated escape occurs when the perpendicular dis- 
tance between the trajectories is larger than the width 
Wl,r of the leads. This requires a minimal time 7V(e)/2 
between encounter and escape, where [3|| 



T w (e) = X-Hnie- 2 (W/L) 2 } 



(32) 



The two pairs of trajectories then escape in an uncor- 
related manner, typically at completely different times, 
with completely different momenta (and possibly through 




(b) Weak-localization and coherent back-scattering for reflection 
correlated r 




FIG. 10: (Color online) Sketches of the trajectory pairing 
which give the leading off-diagonal contributions to (a) trans- 
mission (conductance) and (b) reflection. All contributions 
involve a trajectory 7I crossing itself at an angle e, and a 
trajectory 72 which avoids the crossing. The action differ- 
ence between the two trajectories is thus small and does not 
fluctuate under averaging. For transmission, yo is on L lead 
and y is on R lead, for reflection both t/o and y are on L lead. 
There are two reflection contributions. On the left is weak lo- 
calization, and on the right is coherent backscattering (details 
of the latter are in Fig. I12H . 



different leads). Correlated escape occurs in the other 
situation when the distance between the trajectories at 
the time of escape is less than Wl.r. Then the two 
pairs of trajectories escape together, at the same time 
through the same lead. This latter process affects co- 
herent backscattering (see Fig. I10[). The second new in- 
gredient is that the survival probability for a trajectory 
with an encounter is larger than that of a generic trajec- 
tory. This is so, because the encounter duration affects 
the escape probability only once. In other words, if the 
trajectory did not escape in its first passage through the 
encounter, neither will it during its second passage (this 
was first noticed in Ref. [2?|b 

We calculate the contribution from pairs of transmit- 
ting trajectories sketched in Fig.llOb. For preserved time- 
reversal symmetry, the action difference for this contri- 
bution is SS w i = Epe 2 /X [33. We now note that the 
probability to go from Yo to Y in time t, is the prod- 
uct of the probability to go from Yo to a point on the 
energy surface inside the cavity Ri = (rj , <pi ) (where (j>\ 
defines the direction of the momentum) in time t\ and 
the probability to go from Ri to Y in time t — t-y, when 
one integrates Ri over the energy surface C. Thus the 
quantity P introduced above can be written as 

P(Y,Y ;t) - J dR 2 dR 1 P(Y,R 2 ;i-t 2 ) 

x P(R 2 ,R 1 ;t 2 -ti)P(Ri,Y ;ti). (33) 

where P(R 2 ,Ri;t) is the probability density to go from 
Ri to R 2 in time t, but P(Ri,Yo;i) is a probability 
density multiplied by the injection momentum, pp cos 6q. 
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FIG. 11: (Color online) Sketch of a trajectory (solid line) 
which crosses itself at point ri = (x,y), visiting this point 
first at time ti and second at time £2- Superimposed is an 
infinitcsimally different trajectory (dashed line) which also 
visits the point (x, y) at ti, but is at point r2 = (x + 5x, y + Sy) 
at time ti. This trajectory also intersects itself; however, it 
visits the self-intersection (which is no longer (a;, y)) at times 
ti + Sti and ti + 5t2- 



We then restrict the probabilities inside the integral to 
trajectories which cross themselves at phase-space posi- 
tions R12 with the first (second) visit to the crossing 
occurring at time t\ (i 2 ). Using Fig. ^2 we write dR 2 = 
Vp sinediidt 2 de and set R 2 = (r 2 ,</> 2 ) = {^i,4>i ± e). 
Next, we note that the duration of the loop must ex- 
ceed Ti(e) = A -1 ln[e -2 ], because for shorter times, two 
trajectories leaving an encounter remain close enough to 
each other that their relative dynamics is hyperbolic, and 
the probability of forming a loop is zero. Similarly the 
path cannot transmit unless t\,t — i 2 > IV ( e ); because 
for ti, t — t 2 < ?V(e) the legs (see Fig.|5J) are so close to 
each other that if one leg escapes through a given lead the 
other one will escape with it through the same lead. Thus 
the probability that a trajectory starting at Yo crosses 
itself at an angle ±e and then transmits, multiplied by 
its injection momentum p-p cosQq, is 

roo /•t—T\ V /2 rt 2 —T L 

I(Y ,e) = 2u| sine / dt dt 2 / dii 

JT l +T w JT l +T w /2 JTw/2 

x [ dY [ dRi P(Y,K 2 ;t-t 2 ) 

JR JC 



xP(H^,R i; t2-ti) P(Ri,Y ;ti) 



(34) 



where 7V,Tz, are shorthand for IV (e), Iz,(e). 

To get the weak localization correction to conductance 
we sum only contributions where 7I crosses itself, we 
then take twice the real part of this result to include the 
contributions where 7I avoids the crossing (and hence j2 
crosses itself). Thus 



5wi 



dY deRe[e 



iSS„i/h] 



</(Y ,e)). (35) 



We perform the average of the Ps as follows. Within 
7V(e)/2 of the crossing the two legs of a self- intersecting 
trajectory are so close to each other that their joint es- 
cape probability is the same as for a single trajectory. 



Self-intersecting trajectories thus have an enhanced sur- 
vival probability compared to non-crossing trajectories 
of the same length, i.e. the duration of the crossing 
must be counted only once in the survival probability 
|27). Outside the correlated region, the legs can es- 
cape independently through either lead at anytime. Fur- 
thermore, the probability density for the trajectory go- 
ing to a given point in phase-space is assumed to be 
uniform. Thus the probability density for leg 1 gives 
(P(Ri, Yo;*i)) = (2ttA)~ 1 exp(— ii/r D ) xp F cos9 , and 
the loop's probability density is (P(R 2 ,Rij_£ 2 — = 
(27rA)- 1 exp{-[t 2 -ii-T w (e)/2]/r D } (Ref. gg). Finally 
the conditional probability density for leg 2 (given that 
leg 1 exists for a time t\ > 2V(e)) is (P(Y,R 2 ;i — 1 2 )) — 
[2{W L + Wr)^}- 1 cos 0exp{-[t - ti- T w (e)/2]/r D }. 
Hence one finds 



(P(Y,R 2 ;t-i 2 ) P(R 2 ,R 1 ;t 2 
1 pf cos 9 cos 9o 



(2ttA) 2 2(W L + W k )td 
so that (7(Yo,e)) becomes 



h) P(R 1 ,Y ;t 1 )) 
cxp[-(t-T L (e))/r D ], (36) 



(J(Y ,e)) = — — — Pf sine cos 



nA 



R 



N L + N, 



■exp[-T L (e)/r D ]. (37) 



We insert this into Eq. (|35|) . The e-integral is dominated 
by contributions with e <C 1, so that we write sine ~ e 
and push the upper bound for the e-integration to in- 
finity. The e-inte gral can then be computed to give an 
Euler T- function [23. To leading order in (Atd) -1 it 
equals — 7rfi,(2PFTD) _1 - The integral over Yo yields a 
factor of 2Wl- Finally noting that = (nh)~ 1 pFWi 1 
and (TVl + Nn)^ 1 — (my4.) -1 fiTD, the weak localization 
correction to the conductance reads 



■exp 



l/ru] 



(38) 



We see that weak localization is exponentially suppressed 
with Tg 1 jm in term of the closed cavity Ehrenfest time 



C. Quantum corrections to reflection 

The above result (|38|l has already been derived in 
Ref. [23| by a similar approach. We go beyond that by 
showing explicitly that our semiclassics preserves the uni- 
tarity of the scattering matrix. There are two leading- 
order off-diagonal corrections to reflection. They are 
shown in Fig. HOf b). The first is weak localization while 
the second is coherent backscattering. The former re- 
duces the probability of reflection to arbitrary momen- 
tum, while the latter enhances the probability of reflec- 
tion to the time-reverse of the injection state. The dis- 
tinction between these two contributions is whether or 
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L lead 



cavity 




FIG. 12: (Color online) Trajectories for the backscattering 
contributions to reflection. Trajectory 7I (solid black line) 
start on the cross-section of the L lead at position yo with 
momentum angle #0 and ends at y with momentum angle 6. 
In the basis parallel and perpendicular to 7I at injection, the 
initial position and momentum of path 7I at exit are ro± = 
(yo - y) cos O , r || = (yo - y) sin 9 and p 0± ~ ~p F (0 - ). 



not the legs escape while correlated: for weak localiza- 
tion the legs escape independently, while for coherent 
backscattering the legs escape together in a correlated 
manner. 



1. Weak localization 

The weak localization contribution to reflection, i? w i, 
is derived in the same manner as g w i, replacing however 
a factor of W R /{W L + W R ) with W L /(W L + W R ). One 
obtains 



R 



wl 



Nl 



(JVl + N R )- 



•exp[ 



I/td] 



(39) 



2. Coherent backscattering 

Contributions to coherent backscattering are shown in 
Fig. 110b . with Fig. ^] showing the trajectories in the 
correlated region in more detail. These contributions re- 
quire special care because (i) their action phase difference 
<55 c bs is not given by the Richter-Sieber expression used 
so far and (ii) injection and exit positions and momenta 
are correlated. 

From Fig. El and noting that 726 decays exponentially 
towards 7la, we find the action difference between these 
two path segments to be 



?2fc - Si a = Pf(uo - y) sin d 

+\m\{y - y) 2 cos 2 ( 



'0 ■ 



(40) 



We have dropped cubic terms which only give h- 
corrections to the stationary-phase integral. The action 
difference between 72a and 7I6 has the opposite sign for 
yo — y and 80 replaced by 9. We get for the total action 
difference, in terms of (ro± 1 po±), 



SS C 



-{po± + m\r 0± )r ± 



(41) 



The coherent backscattering contribution to the reflec- 
tion reads 

R chs = {2TTH)- 1 JdY dY J dt 

x (P(Y, Y ; t)) Re [ e i5s ^/ h ] . (42) 

Note that this contains all those contributions where 7I 
crosses itself and all those contributions where it avoids 
crossing (so 72 crosses itself), thus there is no need to 
take twice the real part here (unlike for g w i). To perform 
the average we define T^ v (ro±,po±) and T' L (ro± 1 po±) as 
the times for which the perpendicular distance between 
the 7la and 7I6 is W and L, respectively. For times less 
than T^/(roj_,po±) the escape probability for two trajec- 
tories is the same as for one, while for times longer than 
this the trajectories evolve and escape independently. 
For i? c bs we consider only those trajectories that form 
a closed loop, however they cannot close until the two 
trajectory segments are of order L apart. The i-integrals 
must have a lower cut-off at 2T^(ro±,po±), hence 



/ dY f dt(P(Y,Y ;t)) 

JR J2T' 



Py COS t/o 



N L + N) 



R 



■exp[-T'(r x,Pox)/T D ], (43) 



where T'(r ±,p ±) = 2T' L (r Q ±,p ±) - T^(r ±,_po_L)- For 
small (po_L + rn\ro±) we estimate 



T'(r ±,Po±) ^ A Mn 



W(p ± + m\r ±) 



i\L 2 



(44) 



We substitute the above expression into i? c b s , write 
PFCOS^odYo = dyod(pF sin^o) = dro±dpo± [lj, and 
then make the substitution po = poi_ + m\ro± . We eval- 
uate the ro±— integral over a range of order Wl, take the 
limits on the resulting po-integral to infinity and write it 
in terms of Euler r-functions. Finally we systematically 
drop all terms 0(1) inside logarithms. The result is that 



Rr 



Wl 



N h + N R 



exp[ 



(45) 



Thus we see that coherent backscattering is also sup- 
pressed exponentially in exactly the same manner as 
weak localization. Hence i? c bs + R w \ — — 9wi and uni- 
tarity is preserved. 



D. The off-diagonal nature of coherent 
backscattering 

Continuous families of trajectories that are present in 
open chaotic systems, such as 72 and 7I in Fig.^J have 
an action difference given in Eq. (|41l) . This action dif- 
ference does not fluctuate under energy or sample av- 
eraging, moreover, these contributions are not diagonal 
in the lead mode basis. The stationary phase integral 
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over such trajectories is dominated by po± ~ — m\ro± 
where r ± is integrated over the total lead width. Thus 
Po± varies over a range of order mXW 3> ttH/W, cou- 
pling to many lead modes. Such contributions were not 
taken into account in the previous analysis of coherent 
backscattering [s^. This caused the erroneous belief (of 
the authors of the present article amongst others) that 
coherent backscattering originates from trajectories that 
return to any point in the L lead with 9 ~ ±#o, which 
would have implied that the coherent backscattering was 
independent of the Ehrcnfest time. 

Once we correctly sum the many off-diagonal contribu- 
tions to 



nm Vrira 



2 which have an encounter near the L 
lead, we conclude that coherent backscattering approxi- 
mately doubles the weight of all returning trajectories in 
a strip defined by 



1 - O - -P F 1 m\(y - yo ) cos ( 



(46) 



across the lead. This strip sits on the stable axis of the 
classical dynamics, with a width in the unstable direction 
of order %{p$W) . Therefore, trajectories in the strip 
first converge toward each other, and only start diverging 
at a time of order Tg P /2. Such trajectories cannot form 
a loop on times shorter than r^ p + Tg . 



ti,t 2 gives 



( \2 

(J(Y ,e)) = P\ siiu cusfri. 



irA 



A R exp[-r L (e)/T D ] 
A L + A R 1 + aA 2 (Tu/T { )<i> 2 



(48) 



After a similar analysis for i? c b s , we conclude that for 
finite flux, the quantum corrections to the average con- 
ductance acquire a Lorentzian shape, 



ffwl($) = 

Rwl($) = 

i?cbs($) = 



iV L JVi 



R 



expl-r^/rp] 

(A L + A R ) 2 l + aA 2 (r D /r f )$ 2 

exp[-rg/r D ] 
(A L + N R ) 2 1 + aA 2 (r D /r f )$ 2 



Nr 



exp[-r^/T D ] 



N h +N R l + aA 2 (r D /r f )$ 2 ' 



,(49a) 
,(49b) 
(49c) 



Interestingly enough, there is no Ehrenfest dependence 
in the width of the Lorentzian. 



Weak localization in the two-fluid model 



E. Magnetoconductance 

A weak magnetic field has very little effect on the clas- 
sical dynamics. Its dominant effect is to generate a phase 
difference between two trajectories that go the opposite 
way around a closed loop. This phase difference is Ai oop & 
where Aoop is the directed area enclosed by the loop, 
and <f> is the flux in units of the flux quantum. To in- 
corporate this in the theory we must introduce a factor 
of exp[L4i oop $] into /(Yo,e) in Eq. Q34f) and inside the 
average in f? c b s in Eq. (|4*2~)l . To average (exp[L4i oop <!>]), 
we divide the loop into two parts — the correlated part 
(within Ti,(e)/2 of the crossing), and the uncorrelated 
part (the rest of the loop). We average the two parts 
separately. 

For the uncorrelated part, we use the fact that the 
distribution of area enclosed by classical scattering tra- 
jectories in a chaotic system is Gaussian with zero mean 
and a variance which increases linearly with time |ll| . 
One then has, 



= i^unoo„*\ 



exp[-aA 2 $ 2 (i 2 - h - T L (e))/n], (47) 



where, A is the area of the cavity, a is a system- 
dependent parameter of order unity, and Tf is the time 
of flight between two consecutive bounces at the cavity's 
wall. 

We comment on the correlated part in Appendix [D] 
where we show that it provides at most only small correc- 
tions ©(rf/rp) which we henceforth ignore. Multiplying 
the integrand in Eq. i|34|) with 147( 1 . and integrating over 



Weak localization can also be calculated in the frame- 
work of the special basis constructed in the first half of 
this paper. We can split all contributions to conductance 
into classical and quantum contributions using the classi- 
cal dynamics. By construction the classical modes couple 
to the trajectories shorter than Tg P , while the quantum 
ones couple to the ones longer than Tg P . We cut the time- 
integrals in all the above quantities at Tg P , the result is 
that the classical cavity has 



-^diag 

f7 cl 
i/wl 



At N-a op , 

k [l-e-V/™], 



A L + A R 

N l [l_C-E P /^] ; 



A L + A R 

^wl = ^cbs = 0. 



(50a) 

(50b) 
(50c) 



The quantum modes carrying the remaining contribu- 
tions. Inserting by hand the phase-space splitting into 
the sum rules and bound of time integrations in the 
above semiclassical treatment, one recovers the exponen- 
tial suppression of weak localization, Eq. (|38|l . The quan- 
tum fluid is thus clearly not RMT, which invalidates the 
effective RMT model j4(- This is because contributions 
to weak localization and coherent backscattering come 
from trajectories longer than Tg P + t f , and the propor- 
tion of such trajectories in the quantum cavity goes like 
exp[-T$/TD]. 

The existence of classical bands is key to both the exis- 
tence of the two separate fluids (block diagonal nature of 
S, Eq. (|13|) ~) and to the exponential suppression of coher- 
ent backscattering and the unitarity of the semiclassical 
theory presented here (see also 21]). 



VI. NUMERICAL SIMULATIONS 

We finally check our semiclassical theory for weak lo- 
calization at finite t^/t^, against numerical simulations. 
We consider open systems with fully developed chaotic 
dynamics, for which td > 1. Because r^ 1 grows logarith- 
mically with the Hilbert space size M = , and since 
we want to investigate the regime of finite t^/td, we 
model the electron dynamics by the kicked rotator map 
I2H . The Hamiltonian is given by 



H = 



{p + Paf 



K cos(x + xq) S(t — nTf). (51) 
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The kicking strength K drives the dynamics from inte- 
grable (K = 0) to fully chaotic [K > 7, with Lyapunov 
exponent Arf w \tx{K/2)\. The parameters po and xq 
are introduced to break the Hamiltonian's two symme- 
tries |4lj . Only when these two symmetries are broken 
does one witness a crossover from the j3 = 1 to the 
j3 = 2 universality class corresponding to break- 

ing the time reversal symmetr y lill . The procedure of 
varying po followed in Refs. [2g,|22| results in a strongly 
non-Lorentzian magnetoconductance. Thus the agree- 
ment between numerics and analytics in Ref. [2]} could 
not be extended to the magnetoconductance curve, which 
is the trademark of weak localization. This motivated us 
to perform numerical investigations following the same 
procedure as in Ref. |42j| . i.e. taking a finite, constant po, 
while varying xq. In this case, the magnetoconductance 
curves are Lorentzian. 

We consider a toroidal classical phase-space x,p € 
[0, 27r], and open the system by defining contacts to bal- 
listic leads via two absorbing phase-space strips [xl — 
5x, xl + Sx] and [xr — 5x, xr + 8x], each of them with a 
width 25x = 7t/td- We quantize the map by discretizing 
the momentum coordinates as pi — 2nl/M, 1 = 1,... M. 
A quantum representation of the Hamiltonian l|51f) is pro- 
vided by the unitary M x M Floquet operator U, which 
gives the time evolution for one iteration of the map. 
For our specific choice of the kicked rotator, the Floquet 
operator has matrix elements 



U lv = Me~ im/Mm+lo)2+( - l ' +lo)2] 



(52) 



x e 2mm(l-l')/M e -{iMK/2ir) cos(2;r (m+m )/M) 



with Iq = pqM/2it and tuq = XqM/2-k. 

We restrict ourselves to the symmetric situation with 
^Vr.l = N. A 2N x 2N scattering matrix can be con- 
structed from the Floquet operator U as 01 



S(e) = P[exp(-ie) - U(l - P T P)]~ 1 UP 1 



(53) 



using a 2N x M projection matrix P which describes the 
coupling to the leads. Its matrix elements are given by 



Pn 



1 if n = m S {m\ 
otherwise. 



}UK (i) }, 



(54) 



FIG. 13: Magnetoconductance curves Ap(mo) = g(mo)—g(0) 
for the open kicked rotator model (denned in the text) at fixed 
classical configuration K — 14 (A w 1.95), Td/t{ = 5, and dif- 
ferent Hilbert space sizes, M = 256 (squares, r^ 1 /td — 0.57), 
M = 512 (diamonds, t£/td « 0.64), M = 1024 (upward 
triangles, t^/to « 0.71), M = 2048 (downward triangles, 
Tb/td « 0.78), and M = 4096 (circles, t£/td ~ 0.85). 
The solid line gives the best Lorentzian fit for the M = 512 
curve, Ag{m ) = 0.15 - 0.15/[1 + 1.03 (m /m c ) 2 ]. Data 
have been obtained after averaging over 100000 (1000 clas- 
sically different samples, each with 100 different quasiener- 
gies for M = 256) to 25000 (500 classically different sam- 
ples, each with 50 different quasienergies for M = 4096) dif- 
ferent samples. Inset: rescaled magnetoconductance data, 
A<?(mo)/A(?(2.2m c ), the data collapses onto a single curve, 
confirming out theory that m c does not depend on r^ 1 . 



An ensemble of samples with the same microscopic pa- 
rameters can be defined by varying the position {m 2 ' }, 
i = 1, . . . , N of the contacts to the left and right leads 
for fixed Tp/Vf = M/2N and K [t{ is the time of flight 
through the system, in this particular instance it is the 
time between kicks, see (|51|) ]. We calculate the conduc- 
tance from the scattering matrix, which we numerically 
construct via an iterative procedure as in Refs. 0, H(J . 

In the universal regime r^ 1 /td = 0, Ref. ^3 found, for 
the = 1 to = 2 crossover a magnetoresistance given 
by 



1 



4tt 



1 + (mo/m c 



(55) 



which has to be compared to Eq. (|49a|l . Our task here 
is to investigate the fate of (|55|l as t^/td increases, and 
in particular to check our analytical predictions l|49a|) 
that (i) the magnetoconductance is Lorentzian with (ii) a 
width which is independent of r^ 1 , but (iii) an amplitude 
which is suppressed exponentially oc exp[— t^/td]- 

Fig-EDshows magnetoconductance curves as t^/td in- 
creases while keeping all classical parameters unchanged. 
The data confirm our prediction l|49a|l . i.e. the curves are 
Lorentzian which depend on t^'/td only through their 
amplitude. The inset of Fig. ^5] makes it clear that the 
typical field necessary to break time-reversal symmetry 
does not depend on Tg — after rescaling the magnetocon- 
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M 

FIG. 14: Amplitude of the weak localization correction to 
the conductance as a function of the size of Hilbert space. 
Data have been extracted from the curves shown on Fig. 1131 
either via a fitting of the magnetoconductance curve with the 
Lorentzian Ag(mo) = 7 — 7/[l+7'(mo/m c ) 2 ] (squares) or tak- 
ing 7 = Ag(2.2m c ) (circles). Both methods confirm the ex- 
ponential suppression 7 tx exp[— 7"e /td] (dashed and dotted- 
dashed lines) with an effectively smaller Lyapunov exponent 
A « 1.26 (compared to \n[K/2] » 1.95) 0. 

ductance amplitude, all curves fall on top of each other. 
This was also confirmed by least square Lorentzian fit- 
ting of the magnetoconductance curves with Ag(m ) = 
7 — 7/(1 +7'(mo/m c ) 2 ], which found 7' = 1 ± 0.06 for all 
cases. 

Fig. [21 finally gives a closer look at the suppression 
of the amplitude of the weak localization correction. 
We extracted the amplitude parameter 7 both from a 
least square Lorentzian fitting of the magnetoconduc- 
tance curve, and from the magnetoconductance ampli- 
tude at mo = 2.2m c , g(2.2m c )—g(0). The two procedures 
confirm the exponential suppression of weak localization 
Pfy. in agreement with Ref. 

VII. CONCLUSIONS 

In this paper we investigated the non-coherent trans- 
port properties of open quantum chaotic systems in the 
semiclassical limit. We have shown how to incorpo- 
rate the nonergodic structures appearing in the clas- 
sical phase-space (see Fig. ^) into quantum transport. 
We followed the scattering approach to transport and 
showed how large phase-space structures result in a 
block-diagonal form of the scattering matrix and the 
splitting of the system into two sub-systems put in par- 
allel. One of these sub-systems is of a purely classical 
nature, consisting of deterministic transmission modes 
(transmission eigenvalues are all or 1). We were able 
to calculate the corresponding transmission eigenvectors 
and connect the emergence of determinism to the sup- 
pression of the Fano factor for shot-noise as well as 
the breakdown of universality for sample-to-sample con- 
ductance fluctuations. The classical phase-space struc- 



tures also cause the exponential suppression of coherent 
backscattering, preserving the unitarity of the semiclas- 
sical theory of weak localization presented here. 

At this point, the quantum mechanical subsystem is 
known not to be RMT. Even though all its modes un- 
dergo a certain amount of mixing, and thus carry quan- 
tum effects, weak localization and coherent backscatter- 
ing come from trajectories longer than t^ p + r^ 1 , which 
have an exponentially small relative weight exp[— rg 1 /td] 
in the quantum cavity. The existence of two separated 
fluids is however confirmed. 

We finally point out that the phase-space method de- 
veloped in the first half of this article should work as 
well in regular systems. We however anticipate difficul- 
ties not present in the chaotic systems treated here due to 
the power-law decay of the band areas and diffraction ef- 
fects at the leads. One open question is why open regular 
systems with large dwell times have a RMT transmission 
spectrum l45l 1461 . but a nonuniversal weak localization 
behavior |lltl47|. 

After this work was completed, a preprint appeared 
which reached the same conclusions about coherent 
backscattering |4S| . 
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APPENDIX A: GAUSSIAN WAVEPACKETS AND 
HYPERBOLIC DYNAMICS 

We show that the time-evolution of a Gaussian 
wavepacket in a uniformly hyperbolic infinite system fol- 
lows the Liouvillian flow. This example supports the 
claim that a Gaussian wavepacket in a chaotic system 
follows the Liouvillian flow up to the timescale at which 
the wavepacket becomes so large that the Liouvillian dy- 
namics ceases to be hyperbolic |34| . 

The uniformly hyperbolic Hamiltonian we consider is 
H = p 2 ,/(2m) — m\ 2 x 2 /2. With the change of variable 
p x = (m\/2) 1 / 2 (q + p q ) and x = (2m\)- 1 / 2 (q - p q ), the 
Hamiltonian can be written as 

n = A (qp q - ih/2) . 

Solving the classical Hamilton equations of motion one 
gets q(t) = q(0)e xt and p q (t) = p q (0)e~ xt . Next, it is 
easily checked that a solution of the Schrodinger equation 
is provided by the wavepacket 

(q\m) = ^exp[-A- 2 (t)[ g - q(t) - iA 2 (t) Pq (t)/h} 2 }, 

with A(t) = A(0)e At . Thus we see that a Gaus- 
sian wavepacket remains Gaussian, simply stretched and 
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shifted by the Liouvillian flow. This quantum calculation 
is exact if the system is infinite. In finite systems, initially 
narrow classical distributions undergo a crossover from 
hyperbolic dynamics to diffusive behavior once their ex- 
tension become comparable to some characteristic length 
scale of the system [l|. I n our case we can thus expect 
that Gaussian wavepackets cease to be Gaussian once 
the wavepacket has spread to a width of order the lead 
width. In other words a Gaussian wavepacket will remain 
Gaussian for times shorter than the Ehrcnfcst time |3J| . 



APPENDIX B: ALGORITHMIC 
CONSTRUCTION OF AN ORTHOGONAL 
PHASE-SPACE BASIS 

A basis of coherent states can be made complete (and 
not overcomplete) but not orthogonal, by placing coher- 
ent states at the vertices of a von Neumann lattice • 
This complete basis can be orthogonalized by following 
standard procedures (e.g. Gram-Schmidt orthogonaliza- 
tion), upon which, however, the basis states become very 
different from one another and extended in phase-space. 
We here describe a numerical algorithm which orthogo- 
nalizes a complete basis of coherent states on a von Neu- 
mann lattice, keeping all states identical (up to a trans- 
lation in phase space) while preserving the phase-space 
localization property of those states. 

One starts from coherent states with wavefunctions 

\cs;i,j) = exp [- ||ay| 2 ] exp [aya+jlO). (Bl) 

The creation and annihilation operators are (here h = 1) 

at = 2- 1 / 2 [Q-LP], a = 2- 1 / 2 [Q + LP], (B2) 

where ay = 2~ 1 / 2 [Qi + iPj] and using dimensionless po- 
sition Q and momentum P with [P, Q] = i. The vac- 
uum state, |0), is a Gaussian wavepacket centered at 
P = Q = 0, 

(Q\0) = (Q|cs;0,0) = 7r~ 1/4 exp[-±Q 2 ]. (B3) 

The coherent state IjBlll is a Gaussian wavepacket cen- 
tered at Q — Qi and P = Pj with the same spread in 
both directions, 

(Q\<x;i,j) = Tr-^expIiP.Q-ilQ-Q,] 2 ], 

(P|cs;t,i) = TT-^expHPQ^ilP-P,] 2 ], (B4) 

where we have dropped irrelevant overall phases. 

To get a complete basis, coherent states are placed at 
each vertex but one of a square von Neumann lattice, 
i.e. a regular lattice on the Q-P plane with each unit 
cell covering an area (2tt) |4^. Translational invariance 
means that the empty lattice vertex may be anywhere. 
This basis of coherent states is complete but it is not 
orthogonal. To orthogonalize it, we make the ansatz that 
there exists a set {Pi} such that the wavefunction 

|ps;i,j) = J2 A'/3j'|cs;*' + *. i' + i) (B5) 



obeys 

(ps;fc,Z|ps;i,j) =5ikSji . (B6) 

Note that the form of Eq. (|B5(I is such that we assume 
that the basis-states will be symmetric under interchange 
of Q and P just as the coherent states are. 

To satisfy Eq. (|B6() . we see that the elements of the set 
{Pi} must obey 

Si,o = exp[-(^/2)(z + *'-z") 2 ] (B7) 

i'i" 

We define a set of vectors, {v^}, written in a non- 
orthogonal basis, {e^}, such that — Pi^ a ii, i.e. 

the ith element of the ath basis vector is = 0i- a - 
The basis is chosen such that the basis-vectors have the 
inner product (Si ■ ij) — exp[— (7r/2)(i — j) ]. The con- 
dition that the vectors {v^} form an orthonormal basis 
is 

6 a , = (v {a) ■ v^) = J2 R-M* ■ (B8) 

y 

This is identical to the condition l|B7|) . thus orthogonal- 
izing this set of vectors is equivalent to finding the /3's 
which satisfy Eq. I|B7|) . We use the following algorithm 
to orthogonalize these vectors 

1. Take a complete normalized (but non-orthogonal) 
basis, {vi}. 

2. define a new basis such {v^} such that v\ = — 
\ J2j^i( v i ' v i) v j\- We then choose Ai such that it 
normalizes the vector v[. 

3. Repeat the procedure, taking the new basis {v[} 
and deriving a basis {v'/ }, and so on. 

We take the coherent states described above as the initial 
non-orthogonal basis, so initially = ^ o- In Table Q] 
we present data for the first ten iterations of the algo- 
rithm, by the sixth iteration the results satisfy Eq. I|B7|) 
with an accuracy of < 10~ 7 . Each iteration improves the 
accuracy by more than one order of magnitude. A PS- 
state generated by this procedure is shown in Fig. EI We 
note that Pi decays approximately exponentially with i. 
Thus the PS-states given by Eq. I|B5(I are exponentially 
localized in position and momentum, as shown in Fig. [S] 
Area-preserving stretches, (Q,P,Qj,Pj) — » 
(KQ,K~ 1 P,nQj,K~ 1 Pj), and rotations (Q,P,Qj,Pj) — » 
(Q cos 8 + Psm9,Pcos0 — Qs'm8,QjCos6 + 
Pj sin 9, Pj cos 6 — Q^sin^) are unitary operation 
for any k, 9. Thus the stretched-rotated basis will also 
be orthonormal and complete. This legitimizes the 
procedure discussed in Section IIII Bl for optimizing 
the PS-basis by fitting it to the PS scattering band 
structure. 
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Iterations 
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2 


1.0386424 


-0.1137192 


1.93720E-2 


-3.0781E-3 


5.24E-6 


4.3E-6 


~0 


~0 


~0 


~0 


~0 


4 


1.0357044 


-0.1130793 


1.74448E-2 


-3.0142E-3 


5.478E-4 


-1.024E-4 


1.95E-5 


-3.8E-6 


7E-7 


-IE- 7 


~0 


6 


1.0357044 


-0.1130793 


1.74448E-2 


-3.0142E-3 


5.478E-4 


-1.024E-4 


1.95E-5 


-3.8E-6 


7E-7 


-IE- 7 


~0 


10 


1.0357044 


-0.1130793 


1.74448E-2 


-3.0142E-3 


5.478E-4 


-1.024E-4 


1.95E-5 


-3.8E-6 


7E-7 


-IE- 7 


~0 



TABLE I: Iterations of the orthogonalization algorithm for coherent states on a von Neumann grid. The table shows the 
value of f3i — j3-i as the algorithm is iterated (Pj ~ means that \Pj\ < 10 -7 ). The algorithm converges with accuracy 10~ 7 
after six iterations. Arbitrary large accuracies are obtained with more iterations. The PS-states are found by substituting the 
tabulated values into Eq. I|B5(I . 



APPENDIX C: EDGE-OF-BAND PHASE-SPACE 
STATES 

The tails of the PS-state shown in Fig. decays ex- 
ponentially with the number of lattice points away from 
the center of the PS-state. Strictly speaking, any PS- 
state has thus a finite amplitude outside the band. We 
can however treat PS-states as classical, i.e. completely 
inside one band, if they are more than j max lattice sites 
away from the edge of that band. If we choose j max = 1, 
we would call states "classical" even if they have ~ 3% of 
their squared amplitudes outside the band (this is similar 
to the situation sketched in Fig. 0). If however we take 
jmax — 3, then a PS-state is only classical if less than 
1CP 5 of its squared amplitude is outside the band. The 
number of edge-of-band states (PS-states that are par- 
tially inside, partially outside a band with area > 2Trh c s) 
of a band which exits at time r is 



close to the crossing, i.e. within T^(e/2) of the crossing. 
The situation is depicted in Fig. We consider a loop 
formed after N bounces at the cavity's wall. In the corre- 
lation region, the segment of the trajectory between the 

bounce N-2 
.bounce 3 



A ii where n= 3 



n cob 



(r) ~ 4j n 



W L W K 
2i:h c gL 2 



1/2 



exp[-Ar/2], (CI) 



where K = L, R. Thus the number of edge-of-bands 
states is 



N, 



cob 



K 



Jmax-'* q 



(C2) 



where the sum is over K = L, R and j max is a number of 
order one. Note that the error we make decays exponen- 
tially with j max , hence the choice of acceptable error only 
changes A^ob logarithmically. From Eq. i|C2l) we conclude 
that the edge-of-band PS-states are a subdominant pro- 
portion of the total number of quantum PS-states, and 
can be ignored. 



APPENDIX D: FLUX ENCLOSED BY THE 
CORRELATED PART OF THE LOOP 

Here we analyze the part of the area enclosed by a loop- 
forming trajectory when it is in the correlation region 



bounce 1 
bounce N" 




bounce 2 
bounce N-l 



FIG. 15: A sketch of the area enclosed by the correlated part 
of the weak localization loop. The area An is defined by the 
segment of the loop-forming trajectory between the (n — l)th 
and nth collisions and the segment between the (N — n -+- l)th 
and the (N — n + 2)th collisions, (dashed region). 



(n— l)th and nth collision with the cavity walls, is highly 
correlated with the segment between the (N — n + l)th 
and the (N — n + 2)th collision. We consider the directed 
area A n enclosed by these two segments (dashed region 
in Fig. 115(1 . We assume that \A n \ is uncorrelated with 
|-/4jn#n|) an d take each such area from a Gaussian dis- 
tribution. The typical perpendicular distance between 
the trajectories at time t (measured from the crossing) is 
ztvpeX^ 1 sinh(Anro). Thus we assume {A n } = and 



(A 2 n ) = 



sA" 



(n-l)To 



dt sinh(At) 



(Dl) 



This grows exponentially with n. Thus the sum over the 
An's is dominated by the largest of them with n max = 
Ti(e)/2rf, whose variance is ~ A 2 . Anticorrelations in 
the signs of consecutive directed areas in the correlated 
region further reduce the total directed area. The flux 
enclosed in the correlated region is thus at most ~ A 2 Q 2 . 
This is smaller than the flux enclosed in the uncorrelated 
region by a factor Tf /td <C 1 • 
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